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I.  Findings:  Scientific  Results 

The  PI,  Andrew  Majda,  the  co-PI,  Rafail  Abramov,  and  the  postdoctoral  associates,  Dimitris 
Giannakis,  and  Michal  Branicki  (funded  by  the  DRI)  have  submitted  the  following  papers  all  with  the 
PI  as  author  (and  other  collaborators  listed) 

A)  Mathematical  Techniques  for  Quantifying  Uncertainty  in  Complex  Systems  with  Model  Error  with 
Prototype  Applications 

Development  of  new  uses  of  information  theory  to  quantify  uncertainty,  irreducible  impression, 
sensitivity,  and  long  range  forecasting  skill  (7,8,9,12,  4,6,1  ).  This  work  includes  explicit  simple 
examples  of  irreducible  imprecision  where  imperfect  models  can  be  tuned  to  match  the  climate  mean 
and  variance  of  the  perfect  model;  nevertheless  the  response  to  external  forcing  has  an  intrinsic 
information  barrier  which  cannot  be  improved  within  the  class  of  imperfect  models.  The  development 
and  application  of  these  ideas  to  unambiguous  simple  models  for  turbulent  diffusion  with  complex 
features.  The  theoretical  development  of  algorithms  based  on  fluctuation  dissipation  theorems  for 
sensitivity  and  long-range  forecasting  including  intrinsic  skill  barriers  for  popular  linear  regression 
models.  New  non-Gaussian  filtering  algorithms  for  multi-scale  filtering  of  turbulent  signals 
(14,11,3,2). 

B)  Nonlinear  Laplacian  Spectral  Analysis  (NLSA)  for  Time  Series:  Capturing  intermittency  and  Low 
Frequency  Variability 

Majda  and  Giannakis  have  developed  novel  NLSA  algorithms  and  applied  them  to  comprehensive 
climate  models  (13,10,5).  Many  processes  in  science  and  engineering  develop  multiscale  temporal  and 
spatial  patterns,  with  complex  underlying  dynamics  and  time-dependent  external  forcings.  Because  of 
the  importance  in  understanding  and  predicting  these  phenomena,  extracting  the  salient  modes  of 
variability  empirically  from  incomplete  observations  is  a  problem  of  wide  contemporary  interest.  Here, 
we  present  a  technique  for  analyzing  high-dimensional,  complex  time  series  that  exploits  the 
geometrical  relationships  between  the  observed  data  points  to  recover  features  characteristic  of 


1 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

30  SEP  2012  2' REPORT  TYPE 

3.  DATES  COVERED 

00-00-2012  to  00-00-2012 

4.  TITLE  AND  SUBTITLE 

Report:  Low  Frequency  Predictive  Skill  Despite  Structural  Instability 
and  Model  Error 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

New  York  University, Courant  Institute  of  Mathematical  Sciences, 251 
Mercer  Street, New  York, NY, 10012 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  US 

unclassified  unclassified  unclassified  Report  (SAR) 

4 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


strongly  nonlinear  dynamics  (such  as  intermittency  and  rare  events),  which  are  not  accessible  to 
classical  singular  spectrum  analysis.  The  method  employs  Laplacian  eigenmaps,  evaluated  after 
suitable  time-lagged  embedding,  to  produce  a  reduced  representation  of  the  observed  samples,  where 
standard  tools  of  matrix  algebra  can  be  used  to  perform  truncated  singular-value  decomposition  despite 
the  nonlinear  geometrical  structure  of  the  dataset.  We  illustrate  the  utility  of  the  technique  in  capturing 
intennittent  modes  associated  with  the  Kuroshio  current  in  the  North  Pacific  sector  of  a  general 
circulation  model  and  dimensional  reduction  of  a  low-order  atmospheric  model  featuring  chaotic 
intennittent  regime  transitions,  where  classical  singular  spectrum  analysis  is  already  known  to  fail 
dramatically. 

We  apply  NLSA  to  study  North  Pacific  SST  monthly  data  from  the  CCSM3  and  ECHAM5/MPI-OM 
models.  Without  performing  spatial  coarse  graining  (i.e.,  operating  in  ambient-space  dimensions  up  to 

I. 6  _  105  after  lagged  embedding),  or  seasonal-cycle  subtraction,  the  method  reveals  families  of 
periodic,  low-frequency,  and  intennittent  spatiotemporal  modes.  The  intermittent  modes,  which 
describe  variability  in  the  Western  and  Eastern  boundary  currents,  as  well  as  variability  in  the 
subtropical  gyre  with  year-to-year  reemergence,  are  not  captured  by  SSA,  yet  are  likely  to  have  high 
significance  in  a  predictive  context  and  utility  in  cross-model  comparisons. 

II.  Contributions  to  Discipline 

Development  of  new  uses  of  infonnation  theory  to  quantify  uncertainty,  irreducible  impression, 
sensitivity,  and  long  range  forecasting  skill  (7,8,9,12,  4,6,1  ).  This  work  includes  explicit  simple 
examples  of  irreducible  imprecision  where  imperfect  models  can  be  tuned  to  match  the  climate  mean 
and  variance  of  the  perfect  model;  nevertheless  the  response  to  external  forcing  has  an  intrinsic 
information  barrier  which  cannot  be  improved  within  the  class  of  imperfect  models.  The  development 
and  application  of  these  ideas  to  unambiguous  simple  models  for  turbulent  diffusion  with  complex 
features.  The  theoretical  development  of  algorithms  based  on  fluctuation  dissipation  theorems  for 
sensitivity  and  long-range  forecasting  including  intrinsic  skill  barriers  for  popular  linear  regression 
models.  New  non-Gaussian  filtering  algorithms  for  multi-scale  filtering  of  turbulent  signals 
(14,11,3,2). 

III.  Contributions  to  Outside  Disciplines 

Majda  and  Giannakis  have  developed  novel  NLSA  algorithms  and  applied  them  to  comprehensive 
climate  models  (13,10,5).  Many  processes  in  science  and  engineering  develop  multiscale  temporal  and 
spatial  patterns,  with  complex  underlying  dynamics  and  time-dependent  external  forcings.  Because  of 
the  importance  in  understanding  and  predicting  these  phenomena,  extracting  the  salient  modes  of 
variability  empirically  from  incomplete  observations  is  a  problem  of  wide  contemporary  interest.  Here, 
we  present  a  technique  for  analyzing  high-dimensional,  complex  time  series  that  exploits  the 
geometrical  relationships  between  the  observed  data  points  to  recover  features  characteristic  of 
strongly  nonlinear  dynamics  (such  as  intermittency  and  rare  events),  which  are  not  accessible  to 
classical  singular  spectrum  analysis.  The  method  employs  Laplacian  eigenmaps,  evaluated  after 
suitable  time-lagged  embedding,  to  produce  a  reduced  representation  of  the  observed  samples,  where 
standard  tools  of  matrix  algebra  can  be  used  to  perfonn  truncated  singular-value  decomposition  despite 
the  nonlinear  geometrical  structure  of  the  dataset.  We  illustrate  the  utility  of  the  technique  in  capturing 
intennittent  modes  associated  with  the  Kuroshio  current  in  the  North  Pacific  sector  of  a  general 
circulation  model  and  dimensional  reduction  of  a  low-order  atmospheric  model  featuring  chaotic 
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intermittent  regime  transitions,  where  classical  singular  spectrum  analysis  is  already  known  to  fail 
dramatically. 

We  apply  NLSA  to  study  North  Pacific  SST  monthly  data  from  the  CCSM3  and  ECHAM5/MPI-OM 
models.  Without  performing  spatial  coarse  graining  (i.e.,  operating  in  ambient-space  dimensions  up  to 
1.6  _  105  after  lagged  embedding),  or  seasonal-cycle  subtraction,  the  method  reveals  families  of 
periodic,  low-frequency,  and  intennittent  spatiotemporal  modes.  The  intermittent  modes,  which 
describe  variability  in  the  Western  and  Eastern  boundary  currents,  as  well  as  variability  in  the 
subtropical  gyre  with  year-to-year  reemergence,  are  not  captured  by  SSA,  yet  are  likely  to  have  high 
significance  in  a  predictive  context  and  utility  in  cross-model  comparisons. 
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